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ABSTRACT 

This report explains one method of subdividing 5 
polyhedron into triangular facets and "exploding" 
it into the surface of a sphere. 

A mathematical model is included which explains the 
geometry used in subdividing and transforming the 



TABLE OF CONTENTS 


I. Introduction - - - 

II. Basic Geometry 

A. Definitions 

B. Geometry 

III. Projected Research 

IV. Appendix I ... 

V . Appendi x 1 1 - - - 


1 

3 

4 

24 

25 
48 



i v 



INTRODUCTION 

This report explains one method of subdividing a polyhedron into tri- 
angular facets and "exploding" it onto the surface of a sphere. A struc- 
ture is thereby given which may be used in spherical form. 

The tetrahedron, octahedron or icosahedron are the fundamental geo- 
metrical configurations of the structure. A further subdivision of the con- 
figuration chosen is accomplished by subdividing each principle side of 
each principle polyhedral triangle into any number of segments The order 
of subdivision is determined by subdividing the triangle and the origin (or 
center) of the polyhedron into equal angle segments, using the origin as the 
vertice for subdivision. The points of intersection of the equal angle seg- 
ments with the principle side determines the subdivision along the principle 
side of the principle polyhedral triangle. 

The points of subdivision on each side of the Principle polyhedral tri- 
angle are connected with line segments which are parallel to the two remain- 
ine sides of the principle polyhedral triangle under consideration. They 
intersect at a number of points which define a triangular grid of subdivi- 
sion. Due to the method of subdivision, small triangular "windows" occur 
in the grid. The centers of these windows are found by one of two methods 
(discussed later in this report) and are used as the vertices of a triangular 
grid of subdivision of the principle polyhedral face and are tkr, transformed 
to the surface of the- sphere which circumscribes the polyhedron. The cords 
that connect these transformed vertices thus define the structural grid that 
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comprise the structural configuration desired. 


A mathematical model has been determined which explains the geometry 
used in subdividing and transforming the polyhedron into the structural 
configuration desired. From this model a computer program has been written 
which gives the necessary information needed for construction and analysis 
of the structure. As a further aid in investigation of the various forms, 

d plot routine was developed to give a graphical output of each of the struc- 
tural forms. 


BASIC GEOMETRY 


I 



DEFINITIONS: 

DIHEDRAL ANGLE (f?) an angle formed by two planes meeting in a 
common line. The planes themselves are the faces of the dihedral angle, 
and the common line is the element. To measure the dihedral angle measure 
the angle whose vertex is on the element of the dihedral angle and whose 
sides are perpendicular to the element and lie one in each face of the di- 
hedral angle. 

FACE ANGLE (a) an angle formed by two elements meeting in a com- 
mon point and lying in a plane that is one of the faces of the polyhedron. 

AXIAL ANGLE ( Q) an angle formed by an element and a radius from 
the center of the polyhedron meeting in a common point and the vertex of 
the axial angle sharing a vertex of the polyhedron. 

PRINCIFLE ICOSAHEDRAL TRIANGLE (PIT) any one of the 20 equal 
equilaterial triangles which forms the face of the regular icosahedron. 

PRINCIPLE SIDE any one of the three sides of the .principle icosa- 
hedral triangle. 

FREQUENCY the number of equal parts into which a principle side 
is subdivided. 

FACES the triangles making up the "exploded" structural r'orm. 
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GEOMETRY 

Tnis program works with a tetrahedron, octahedron, or icosahedron 
circumscribed by a unit sphere. The icosahedron was chosen as an example 
to illustrate the geometry of the program. The icosahedron is oriented 
in a three dimensional rectangular coordinate system so that the vertices 
of one principle icosahedral triangle are: 


(x , y , z ) 

= 

(o, s/T , 

1 ) 



yr 

{/b/T 


a 

(0, .850651, 

.525731 ) 

(x , y , z ) 

= 

(1 , o 

* v / T ) 



s/5s/t 

yr 



(.525731. 0. 

.850651) 

(x , y , z ) 

= 

( /T , 

1 , 0) 


yr ys/r 

= (850651, .525731, 0) 

where = 1 + yT 

2 

with the intersections of the axis X, Y, Z, located at the origin (0,0,0) 
of the icosahedron, Figure 1. 
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Fiqure 1 

This principle icosahedral triangle is divided into smaller triangular 
units which are "exploded" onto the surface of a sphere constituting 
the desirable space form. 

Using the following formula the planes consisting of the edges of 
the PIT and the origin (Xj, Y lt Zj (X 2 , Y 2 , Z 2 ) (X 3 , Y 3 , Z 3 ) are 
rotated from 3-space into 2-space, Figure 2, 
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Moure 2 


x = AjX + y,y + v 2 z 

y" = a 2 x + u,y + v 2 z [i] 

z' = A 3 x + y 3 y + v 3 z 


Where A, y, v are direction cosines of the X'-axis, Y'-axis, and Z'-axis 
respectively with respect to the old axis and are found by: 

2 7 

Al = x i / / x i + yi + Zi 




^2* ^ 3 * U 2 • 


u 3 ; and v 2) v 3 are found similarly. 


r 


* 


The edge of the PIT is subdivided into units by the following method. 
Figures 3 and 4. 


FIND: the angle 4) contained within the rotated triangle consisting 

of PjPa, and the origin with the vertex located at the origin. 


4 > 

where r 


Arc tan ( 



r 


P x> 

1 and is considered constant 


[ 2 ] 


Y 



THEN: subdivide the angle <p into N angles e 

6 ’ _±_ . T 
N 

where T = Increment 1 to N 


Y 



[3] 
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The points of intersection of 0F 3 and P 1 P 2 are found: 


P : P 2 is 


0P 3 is 


y - Yi 


V 2 - Vi 


x - x, 


y - 0 


x 2 - x t 


y 3 - 0 


[ 4 ] 


x - 0 


x 3 - 0 


The equation takes the following form: 


P1P2 is x(y 2 - yj + y(xi - x 2 ) = yi(xj - x 2 ) + 

Xj(y 2 - yi) 

let (v 2 - yi ) = aj 

(x 2 - Xi) = bi 

Yj (Xj - x 2 ) + x t (y 2 - yj = Cj 


OP is xy 3 - yx 3 = 0 


let 


y 3 = a 2 

-x 3 = b 2 

0 = c 2 


Solve the equations for the point of intersection: 



a i 



Cl 

bi 

D = 



X = 




*2 

b a 


C 2 

b 2 


r w 


Y B Bl Cl 
a 2 c 2 

D 

Rotate the points of intersection along the PIT edge from 2-spaces back 
to 3-spaces. 

x = XjX" + y x y' + VjZ 

y = A 2 x' + \i 2 y" + v 2 z [ 6 ] 

z = A 3 x' + y 3 y' + v 3 z 

where A, y, v, are direction cosines of the X"-axis, Y'-axis, and Z'-axis 
with respect to the old axis and are found: 

A = x / /x 2 + y 2 + z 2 

1 i v i J 1 i 

u i s *i //*7 + v * v 

v . * *. v *, 2 * V * V 

^ > X ; y , y ; and v , v are found similarly. 

2 3 2 3 2 3 J 

Retain the co-ordinates along the edges S , S and S as shown in Fiqure 5. 

12 3 3 
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Figure 5 

After finding the unit measurements along the edges of the PIT, they 
are connected through a 3-way grid determining a smaller triangular grid 
network. Since the units along the PIT edge are not equal, the 3-way grid- 
ding will create "windows' . The centers of these "windows" must be found to 
establish the final 3-way grid network, Figure 6. 


S, 

P 2 ' Pu 



Figure 6 



The gridding and windows are found by the following method: From the 

coordinates along the edges of Sj,S 2 and $ 3 , calculate coordinates of the 
window by finding the intersection of P 3 P 2 with P 3 P 4 and P7“P 2 with P7~P 6 
and Py P 4 with P' 5 P 6 by using the two point form of the equation of a line 
in three-space for the three lines and solve simultaneously for the points 
of intersection. 


x - Xj y ■ y 3 z - 2 1 


v x i y,-y! 


V z i 


P 3 P 4 is: 


X - x 3 y - y 3 . z - 2 


[7] 


x .' x 3 y4-y 3 -4 -a 


Z, -Z. 


P 5 P 6 is: X - x y - y 5 Z - z, 


X 6 ” X 5 


y 6 -y 5 


Z 6- Z 5 


To find the intersection of r 3 P 2 with P 3 P 4 the equation takes the following 
form: 

(1) is; x(y 2 -yj + y(x r x 2 ) = y 1 (x 1 -x 2 ) + x 1 (y 2 -y 1 ) 

(2) P^P, is: y(z 2 -z 1 ) + z(y 2 -y 3 ) = z 1 (y 1 -y 2 ) + y 1 (z 2 -/ 1 ) 

(3) FJP 4 is: x(y 4 -y 3 ) + y(x 3 -x 4 ) = y 3 (x 3 -x 4 ) + x,(y 4 -y 3 ) 

(4) F^ 4 is: y(z 4 -z 3 ) + z(y 3 -y 4 ) = z 3 (y -y 4 ) + y (z 4 -z 3 ) 
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For P ; P 2 let: (y 2 -yj = a, 

(x r x 2 ) = b j 

yJXj-Xj) + x,(y 2 -y,) = c t 

For Pp 4 let: (y 4 -y 3 ) = a 2 

(x 3 -xj = b 2 

y,(x 3 -x 4 ) + Xj (y 4 -y 3 ) « c 2 

using the formula [6] solve for x and y coordinates of the Intersections of 
P? 2 with P^ 4 . 

Find the z coordinate: 

For P^ 2 let: (z 2 -z,) = a t 

(yry 2 ) = bj 

z i(yi-y : ) + yiU 2 -z 3 ) = c 3 

For P 3 R let: (z 4 -z 3 ) = a 2 

(y 3 -y 4 ) “ b 2 

z 3 (y 3 -y 4 ) + y 3 (z 4 -z 3 ) = c 2 

The other two vertices of the window are found in a similar manner. 

Once the coordinates for the vertices of the window are determined, its 
center is found by one of the following two methods: 

METHOD I: 

On the Principle Icosahedral Triangular Plane the 
windows appear as equilateral triangles with vertices 
P ; (x 3 y i* j ) , P 2 (x 2 y 2 z 2 ), P 3 (x 3 y 3 z 3 ) as shown in Figure 7. 
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Figure 7 


METHOD II: 


The center C(cx,cy,cz) is found with the following 
formul a : 

CX = X]+x 2 + X3 
3 

CY = y ' +y ^ +y 3 [8 ] 

3 

CZ = ZJ+Z.+Z3 

3 


The coordinates of the window found on the surface 
of the PIT are first "exploded" to the surface of the 
sphere. The center of the exploded window is then 

L3- 



found by the intersection of angle bisectors. 

To find the projection of each vertex of the window 
onto the unit sphere, translate each vertex along a 
line through the vertex of the PIT and the origin; 
each coordinate of each vertex, PIT, is divided by the 
distance between the vertex PIT and the origin, 

Figure 8. 


y 




[ 11 ] 


P1 TX 

■ Pix 

- P 3 X 

P t -ry 

* p iy 

- P 3 y 

P 1TZ 

• 

- P ,Z 

P 

2TX 

P 2X 

P 3 X 

p 

2 T y 

= P 2X 

- P 3 X 

p 

' 2TZ 

= P 2Z 

' P 3 Z 

p 

3TX X 

P 

’ 3Ty,» 

P 3TZ 


Rotate plane P l5 P 2 , P 3 so that P,P 3 will fall on the X-axis and P 3 is at 
the origin using equation [1]. 

The center is found with the intersection of two angle bisectors of the 
triangular window P,P 2 P 3 , Figure 10. 


Y 
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The angles y and 5 are found: 


Arctan y 2 

— Y 

X 2 

Arctan y, 

2 =, 6 

x r x 2 

rotate P 2 about P 3 toward P 5 , 1$ y degrees 

x 4 = x 2 COS v + y 2 sin Sj Y 

y 4 = y 2 cos >5 y - x 2 sin ^ y 

locate the origin at P lt then rotate P 2 about Pj toward P 3 ig 6 degrees. 
x 5 = (x 2 -Xj) cos ij 5 - y 2 sin % 6+ x 1 
y 5 ~ 'Jz cos fi+ (x 2 - Xj) sin 6 

thus defining pTp 5 and FTP 4 . 

With P 3 at the origin formula l 4 ] may be used to solve for the inter- 
section of line P 3 P 4> PjPj finding center C. Rotate C back to three space 

using formula [6]. Then translate center C back to three space ("C 1 is 
located in the previously 'exploded" window). Figure 11. 

C"x = Cx + P 3 x 

C y = Cy + P 3 y [15] 

C"z = Cz + P 3 z 


[ 12 ] 


[ 13 ] 
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For Method I or Method II, the centers found are "exploded" to the 
surface of the sphere using formula [9] and formula [10]. 

Using the coordinates, this program finds the lengths of the elements 
of the structure (£•), the angle between pairs of elements (face angle a), 
the angle between the elements and a radius from the origin to an endpoint' 
of the element (axial angle ft), and the angle between adjacent faces of 
the structure (dihedral angle 6), Figure 12. 
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Figure 12 


To find the angle between elements the face i a, we use the coordinates 
of their endpoints. The vertex of the angle is a common endpoint to each 
element and is translated to the origin. The other two endpoints Pj and P 2 
are translated in the same manner. Letting (x b y u Zl ) and (x 2 , y 2 , z 2 ) 
be the points resulting from the translations of the endpoints Pi and P 2 , 


cos a 

where dj 
and d. 


x ! x 2 + y^a + z : z ; 


d,d 

y x i 2 + yi + z i 


l i. 

2 2 T 


2 2 2 
x 2 + y 2 + z 2 


a is the desired angle. 

To find axial angles the above method is used except that the vertex 
is established at one end of an element and the origin is used with the 
other endpoint to define the angle. The desired angle is 52. 

The angle between two adjacent faces, the dihedral ? 0, is found 


using 


cos 8 


- AjA 2 + BjB 2 + CjC 2 

A' + b, ! * c,Va 2 2 ♦ Bi ! 



where 


6 is the desired angle. 

A i X + BjY + C,Z + Dj = 0 defines trie plane cortaining one face and 
A 2 X + B 2 Y + C 2 Z + D 2 = 0 defines the plane containing the other face. 
The negative sign is used because the obtuse angle is desired. 

The A, B, and C for each plane are computed as 



1 


x i Y, 1 

c = X 2 Y 2 1 

x 3 Y s 1 

where (X t , Y t , Z 3 ), (X 2 , Y 2 , Z 2 ), and (X,, Y 3 , Z 3 ) lie in the plane. In 

particular the three vertices of each face are used. 

The length of the elements £ are found by using the general equation: 

i = v/Tp - P ) 2 + (p ~P ) 2 + (P ~P ) 2 
v x i x 2 >'i y 2 z i z 2 

£ is the desired length. 

To reduce total output, this program takes into account certain sym- 
metries and outputs only a part of the total angles and lengths. The rest 
of the values are the same as at least one outputed value and can easily be 
found using the following symmetries, Figure 13. 
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FACE ANGLES 


For every face angle opening directly towards (or away from) the point 
(0,0), there are equal angles opening towards (or away from) the point (N,0) 
and (N,N,. For example, the angle (1,1), (0,0), (1,0) with vertex at (0,0) 
is equal to the angle (N-1,0), (N,0), (N,l) and the angle (N,N-1), (N,N), 

(N-l, N-l). Thus, only the face angles facing directly towards or away from 
(0,0) are computed. If the vertex is to lie at (I,J), the angle will be 
either (I + 1, J + 1), (I, J), (I + 1, J) or (1-1, J-l), (I,J), (1-1, J). 

Also, only the face angles falling on the right of or on a line passing through 
(Xj, Yj, Zj,) and the midpoint of the opposite side are computed. 

The elements of the structure can be put into one-to-one correspondence 
with the lengths and dihedral angles. The dihedral angle associated with an 
element is the angle between the two faces containing the element. For each 
element, there are two axial angles, one at each end, but since the element 
is a cord of the circle, the two angles are equal and may be considered one. 

In this case, we have a one-to-cne correspondence between elements and axial 
angles. This program will only compute values around elements parallel to the 
side opposite ( X x , Y x , Zj and on the right side of a line through (X t , Y t , 

Zj) and the midpoint of the opposite side. All other lengths and angles are 
symetric to one of the lengths and angles computed in this manner. 
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PROJECTED RESEARCH 

Another method is being developed of subdividing a polyhedron and 
"Exploding" it into a spherical forir. A correlational study will be con- 
ducted to determine the difference in material usage, volume, and surface 
area of the subdivided forms, using the method described in this report, 
previously developed methods, methods presently under development, and 
the Fuller methods of subdividing a sphere. 

A mathmatical model is also being developed which will explain the 
transformation of polyhedral forms as expandable structures. 
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APPENDIX I 

The computer program here contained, was written for the IBM 7044 
computer, utilizing FORTRAN language. The program may be used for a 
Tetrahedron, Octahedron, or Icosahedron, depending upon the coordinates 
used as input data. The output is given in units based upon a radius of 
1 for the spherical form and therefore may be used as a basis for deter- 
mining large structures. 

The example of the output data given here is for a six frequency 
icosahedral sphere, and may be read as example I. The output takes advan- 
tage of symmetries within the spherical icosahedron as discussed in the 
text material. See Figure 14. 
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CHART 61 


61 Read the parameters for a DO loop governing the frequencies and two 
variables; one indicating the basic structure (tetrahedron, octa- 
hedron, or icosahedron) and another indicating what method of division 
will be used. 

C1-C2 If the card was blank (first parameter zero), s top. 

D1 Determine from the card read, what coordinates are to be used as the 
basic structureand the coefficients of some equations. This informa- 
tion is stored in an array by a DATA statement. 

El Enter a DO loop governing the frequency. 

FI Set the summations for total area, length, and volume to zero. There 
are two summations of length: one for total structure and one above 
a face of the basic structure. 

61 Write a leader record containing the name of the basic structure, the 
frequency, and the method of division. 

HI - J1 Find and write the total number of faces, edges, and vertices both 
in the entire structure and above one face of the basic structure. 

These figures are used for display only and are not used in any j: urther 
calculations. 

K1 Enter DO loop to determine the first end point of the edge. 

A3 Enter DO loop to determine the second end point of the edge. 

B3 The end points of the edges should not be the same. 

C3 Rotate the plane defined by the end points and the origin into the X-Y 
plane with one of the end points into the X axis, and calculate the 
angle <t> = ARCTaN (y'Vx 2 ) necessary for subdivision. 

D3 Enter DO loop to determine the subdivision of the edge. 

E3 6enerate points of subdivision along the edge by N subdivisions of angle 

4 > • 

F3 Rotate the generated points back into the original plane. 

63 Store the points generated into their respective arrangement. 

H3-K3-J3 If DO loop is complete there should be no more points on their 
respective edges. 

A4 The three vertices of the oasic structure are arbitrarily ordered. The 
first vertice is exploded to the surface of the unit sphere to explode a 
point. This is done by dividing the coordinates of a point by its dis- 
tance from the origin. 
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B4 Print out the coordinate of this point: G2J1 -G2K1 ; 62E2-G2F2; 

G2C3-G203. 

C4 G1J5-G1K5; G2F1-G2G1; these blocks are used to count lines printed 

and to provide margins at the top and bottom of each page. 

D4 Subdivision of the basic face will be generated by using the coordinates 

stored in SD1 , SD2, and SD3 and taking all lines parallel to their 
respective edges. The intersections of these lines creates windows 
of which their coordinates are necessary for finding their centers (two 
meth^us of finding the windows center are used) for completion of the 
subdivision. Enter DO loop which generates points along the edge from 
the first vertex to the second. 

E3 Enter a DO loop which generates points along the edge from the first 
vertex to the second. These points will be the first in each row of 
points parallel to the edge from the second vertex to the third. The 
first row will be the one nearest the first vertex. 

E4 Enter DO loop which generates the points alone the Tow determined in 
E4. 

F4 Determine the I. 2. values of the basic unit. 

G4 Enter DO loop which generates the coordinates of the basic unit. 

H4 Set up constants of points in basic units used in determining coordinates 
of the windows. 

J4-K4 All points should be within the PIT. If not, set equal to zero 
(necessary for finding dihedral angles along PIT edges.) 

A5-B5 Points along the PIT edges may be found directly from an SD array. 

C5 INTERC i s a subroutine that finds the coordinates of the intersection 
of 2 lines; used to find the coordinates of the vertice of the window. 

D5 Determine the center of the window by one of two methods; Method #1: 

find the center by taking the intersection of four lines, each passing 
thru a vertice of the window and the midpoints of the opposite side. 

E5 This is the end of the DO loop that finds all of the points within the 
basic unit. 

F5 This explodes all points of the basic unit to the surface of the sphere. 

H5 Writes the coordinates of Pj and P. 
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CHART GII 


B1 


D1 

El 

HI 

J1 

B2 

C2 

D2 

E2 

H2-K2 

A3 


B3 

E3 


G3 


H3-J3 


Due to symmetry, not all angles and lengths are needed. 

Find the length between P l and P 2 . 

Add a multiple of this length into the sums. A multiple is used 
because other elements symmetric to it will not be found due to the 
conditions of B1 . 


Lable and write this length. 

Find the axial angle with vertex angle at P. and legs through P, and the 
origin. 

Lable and write this axial angle. 

Find the total area of and volume under the face generated by P,, P,. 
and P 3 . J 1 2 

Add this into the summations. 


Find the face angles with vertices at P 2 and P 4 and legs passing through 
Pj and P 3 , 

Lable and write the face angles. 

If Pi is in the last row of the face of the basic structure, then 
there is not face between P lf P 3 , and P 4 . Otherwise, the area of and 
volume under that face will be added into the summations. 

Find the dihedral angle between the two faces formed by P 3 , P 2 , p 3 , and 
P„. If P 4 is not above this face of the basic structure, it will not 
be properly generated. In this case, find the angle between the face 
formed by P x , P 2 , and P 3 , and a plane passing through P 3 , P 3 , and the 
origin, doubling this result. 

Lable and write this dihedral angle. 

This is the end of the DO loop of block GIE4. If the loop hasn't 
been satisfied, the next point will be the first in the next row. 

This is the end of DO loop of block GID4. If the loop hasn't been 
satisfied, the next point will be the first in the next row. 

Find, lable, and write the total areas, volumes, and lengths of the 
structure. 


K3 


This is the end of the CO loop of block GIE1. If the loop hasn't been 
satisfied, the same basic structure will be used for the next frequency. 
Otherwise a new card is read. 




1 


As used in blocks GIF! , G1H1-G1J1: 

the vertices above a face are those vertices exploded from points on 
that face. The lengths and areas above the face are those created by 
the vertices above the face. 


As used in blocks G1H1-G1J1, G1E5: 

the basic unit consists of four points combined in two adjacent tri- 
angles with two parts approximately parallel to the edge containing 
vertex 2 and 3. 
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MBFTC GEOD NODECK 

DI MENS I CN SD 1(3, 32), SO 2(3, 32), SD 3(3, 32) • P|4, 31 
DIMENSION COORD ( 3, 3, 3 ), V(3, 3), LAB ( 2 , A), NAMES! 2, 3), 

1 F I G ( 2 , 3) 

INTEGER EFV ( 3 » 3), PAG SiZ, TE OC IC, VL , EL, FL, VG , EG, FG 

REAL L OF TRI, LENGTH, L CF FIG 

DATA FIG / 4», 8*, 6., 16*, 20 »» 40. / 

DATA NAMES / 12H TETRAHEDRON, 12H OCTAHEDRON , 12K ICOSAHEDRON /« 
1 PAG SIZ / 56 / 

DATA CCORO / -.57735C27, -.57735027, .57735027, .57735027, 

1 -.57735027, -.57735027, -.57735027, .57735027, -.57735027, 

2 1., 0., 0., 0., 0., 1., 0., 1., 0., .0, .P5C65Q81 , .52573111, 

3 .85065081, .52573111, 0., .52573111, 0., .85065081 / 


DATA EFV / 

4 , 8 , 

20, 2, 

4, 

10, 6, 12, 

30 / 

EQUIVALENCE 
1 (V(l, 1), 

XI ), 

(V (2 , 

1), 

Yl), C V I 3 , 

1), Zl), 

2 (V(l, ?), 

X 2 ) , 

( V ( 2 , 

2), 

Y 2), ( V( 3, 

2), Z 2 ) , 

3 (Vtl, 3), 

X3) , 

( V ( 2 , 

3), 

Y3 ) , ( V( 3, 

3), Z3), 


4 

( LAB ( 1 , 

1 ) 

, ID, 

( L AB( 2 , 

D, 

Jl), 

(LAB ( 1 , 2), 

12), 

5 

( LAB (2 , 

2) 

, J2 ) , 

(LABI 1, 

3), 

13), 

(LAB (2, 3), 

J3) * 

6 

(LABI 1, 

4) 

♦ 14), 

( LABI 2, 

4), 

J4) 



COEF ( A 1 , 

A2, 

A3, B 1 

, E2, B3 ) a 

A1 * 

( B2 - B 3 ) - 

A2 * ( B 1 - B3) 4 


1 A3 * (81 - 82) 

70 FORMAT (A 12) 

71 FORMAT ( 10HIFREQUENCY 15, 2X, 2A6, 5X, 8HMETH0D 2) 

72 FORMAT (7H0V(L) = 15, 9H E(L ) » 15, 9H F ( L ) « 15, 9H V(G) - 

* 15, 9H E(G) « 15, 9H FIG) » 15 ///) 

73 FORMAT (9HL 0 0 3(5X, F10.6)///) 

74 FORMAT ( 1HL 214, 3(5X, FIG. 6)) 

75 FORMAT £1HA) 

76 FORMAT (9X, 8HLENGTH , 2I3X, 214), 21X, F11.8) 

77 FORMAT (9X, 8HAXIAL , 11H 0.0 0.0, 2(3X, 214), 10X, 2F11. 6) 

78 FORMAT (9X, 8HFACE , 3( 3X, 214), 10X, 2F11.6) 

79 FORMAT (9X, 8HDIHEDRAL, 2(3X, 214), 21X, 2F11.6) 

80 FORMAT (//) 

81 FORMAT (12H1 CNE FACE/ 9X, 4HAREA, 27X, F15.6, / 9X, 6HV0LUME, 

* 25X, F15.6 / 9X , 6FLENGTH, 25X, F15.6 // 22H COMPLETE STRUCTU 
*RE / 9 X , 4HAREA, 27X, F15.6 / 9X, 6HV0LUME , 25X, F15.6 / 9X, 

* 6HLENGTH, 25X, F15.6) 

1 READ (5, 70) MIN, MAX, INC, TE OC IC 
IF (MIN .EQ. 0) STOP 

IF (MAX .EQ. 0) MAX » MIN 

IF ( INC .EQ. 0) INC » 1 

K * 1 

IF (TE OC IC .EC. 8) K * 2 
IF (TE OC IC .EQ. 2C ) K 3 3 
DO 2 I * 1, 3 
DO 2 J * 1, 3 

2 V(J, I) * COORD ( J , 1, K) 

DO 25 N * MIN, MAX, INC 

FLOAT N * N 

SUM OF L a o. 

SUM OF A a o. 

SUM OF V * 0. 

L OF TRI * 0. 
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WRITE (6, 71) N, NAPES! 1 , K), NAMES ( 2, K) 

VL = ( N ♦ 1 ) * ( N ♦ 2 ) / 2 

EL = 3 * N * IN ♦ 1) / 2 

FL = N**2 

VG = EFV ( K t 2) * N**2 + 2 

EG = EFV IK, 3) * N**2 

FG = EFV (K, 1 ) * FL 

WRITE (6, 72) VL, EL, FL, VG, EG, FG 

LI = 0 

DO 5 L2 = 1, 2 
DO 5 LN = 2, 3 
IF (L2 .EQ. LN) GO TC 5 
LI = LI + 1 

CALL ROTATE (l, VII, LN), V<2, LN), V{ 3, LN )» X2R , V2R, Z2R, 
1 VII, L2 ) , V ( 2, L2>, V( 3, L2)) 

THETA = ATAN2(Y2R, X2R ) 

DO 5 14 * 1, N 
T = L4 

ANG = T * THETA / FLOAT N 
X = COS (ANG) 

Y = SIN (ANG) 

A1 = -Y 

81 = X 

A2 = -Y2R 

82 = X2R - l. 

C2 = 62 

Y = ( A1*C2 ) / ( A 1*B2 - A2*B1) 

X = (B1*C2) / ( A2*B l - A 1*82 ) 

CALL ROTATE (2, X, Y, 0.,XR, YR, ZR ) 



IF 

(LI 

- 2) 

27, 

27 

SD 

1(1, 

L4) 

= 

XR 


SD 

1(2, 

L4 ) 

= 

YR 


SD 

1(3, 

L4 ) 


ZR 


GO 

TO 5 




3 

SD 

2(1, 

L4 ) 

= 

XR 


SD 

2(2, 

L4 ) 


YR 


SD 

2(3, 

L4 ) 

= 

ZR 


GO 

TO 5 




4 

SD 

3(1, 

L4 I 

= 

XR 


SD 

3(2, 

L4 ) 

S 

YR 


SD 

3(3, 

L4 ) 


ZR 


5 CONTINUE 

R = SQRT ( X 1**2 ♦ Y 1 **2 ♦ Zl**2) 

P(l» 1) = XI / R 
P( 1, 2 ) = Y 1 / R 
P( 1. 3) = Zl / R 

WRITE (6, 73) P(i, 1), P| I, 2), P(l, 3) 
LINES = 12 
DO 24 II - 1 , N 
DO 23 J1 = l, II 

12 = II - 1 

J 2 * Jl - 1 

13 * 11 

J 3 * J2 

1 4 * II + 1 
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A 


J4 

X 

J1 


00 

11 

LA = 

1. 4 

LI 

X 

LAB( 1 , 

LA) 

LJ 

= 

LAB (2 , 

LA) 

IJ 

= 

Ll - LJ 


NJ 

X 

N - LJ 


IF 

(LI .LE. 

N) GO TC fc 5 

P( LA, 

1) = 0 

• 

Pf L A , 

2) = 0 

• 

P ( L A , 

3) = 0 

• 

GO 

TO 

11 



65 IF (LI .EQ. LJi GC TC 7 
IP (LJ .EQ. 0) GO TC 9 
IF (LI .NE. N) GO TC 10 
PI LA, 1 ) * SD 3(1, NJ ) 
P(U, 2) = SO 3(2, NJ ) 
P(LA, 3) = SO 3(3, NJ ' 
GO TO 11 

7 IF (LJ .EQ. 0) GO TC 8 
P( LA, 1 ) = SO 1(1, LI ) 


SO 2(1, LI ) 
SO 2(2, LI) 
SO 2(3, LI ) 


P ( LA , 2 ) = SD 1(2, LI) 

P (LA, 3 ) = SD 1(3, LI) 

GO TO ll 

8 P ( LA , 1) = XI 
P (LA, 2) = Y 1 
P ( LA, 3) = Z1 
GO TO 11 

9 P ( LA , 1 ) = SO 2 ( 1, LI ) 

P(LA, 2) = SO 2(2, LI) 

P(LA, 3) = SO 2(3, LI ) 

GO TO 11 

10 CALL INTERC (SO 1(1, LI), SO 1(2, LI), SD 1(3. LI) 9021. ... 

\ cn i <2# LI>? S ° 2<3f LI,# S0 3flt SD 3C 2 1 N J) • SD *3 ( 3 

2 SO 1(1, LJf, SO 1(2, LJ), SO Us! LJ), XEl. YEl.2lh) ’ * 

i^^cn *^D 3(1, IJ), SO 3(2, IJ), SO 3(3, IJ), SO 2(1 IJ ) 

p cn f ! 2 ’ SD 3(1, NJ), SO 3(2, NJ) t SO 3(3 HJ) 

2 SD 1(1, LJ), SO 1(2, LJ), SO 1(3, LJ>I XE2, YE2. ZE2) * 

A <n ,^ ERC <S ° 3 * 1 ’ * J, » SD 3<2 » IJ), SD 3(3, IJ), SO 2(1, IJ), 

1 SD 2(2, IJ), SO 2(3, IJ), SO 1(1, LI), SO 1(2. LI) SO i u itt 

2 SD 2(1, LI), SO 2(2, LI), SD 2(3 L XE3, VEst 1 e 3 ? * 

?<ti: ?r! R r^ri c ? e r. c i^r/°L w,NDow ,n F .« . . 

P 1 LA, 2) = ( YE 1 ♦ YE2 ♦ YE3) / 3. 

P(LA, 3) = (ZE1 + ZE2 ♦ ZE3 ) / 3. 

11 CONTINUE 
L = 4 

IF ( 14 .GT. N ) L = 3 
00 12 INDEX 1 * 1, L 

So^^ImJT^uV'** 2 * •’ ,IN0EX »• 2 >**2 ‘ P ( INOEX 1. 

12 P ( I NOE X 1, INDEX 2) » P< INDEX 1, INOEX 2) / R 
IF ( J 3 .NE. 0? GO TC 125 

lines J3 ’ ,,,3 • u * p,i ’ 2I> p ' 3< 31 

IF ( L I f v E S .LT. PAG SIZ) GC TO 125 
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WRITEIo, 75) 

LINES = I 

125 WRITE (6, Z4) II, Jl, P(i, 1), P | 1 , 2), p ( 1 , 3) 

LINES = LINES 4 3 

IF (LINES .LT. PAG SIZ) G r TO 13 

WRITE (6, 75) 

LiNES = 1 

13 IF (2 * J 3 .GE. 13) CC TC 23 

LENGTH = S QRT ( ( P ( 1 , 1) - F ( 3 , I ) ) **2 4 (P(l, 2) - P(3, 2»)**2 ♦ 

* ( P( I , 3 ) - PC 3, 3) )**2) 

L CF TPf = L GF TRI ♦ 3. * LENGTH 

IF (2 * JI .GT . ID GC TC 14 

L OF TK I = L OF TRI 4 3.* LENGTH 

SUM CF L = SUM CF L t 1.5 * LENGTH 

IF (II .LT. N) SUM CF L = SUM OF L + 1.5 * LENGTH 

GO TJ 15 

14 SUM OF L * SUM OF L 4 .75 * LENGTH 

IF (II .LT. N) SUM CF L = SUM OF L ♦ .75 * LENGTH 

15 WRITE (6, 76) 11, Jl, 13, J3, LENGTH 

LINES = LINES + l 

IF (LINES .LT. PAG SIZ) GC TO 16 


WRITE (6, 

75 ) 



LINES = l 




16 XT = P( 1, 

1 ) - 

P( 3, 

1 ) 

YT = PI 1, 

2) - 

P<3, 

2 ) 

ZT = F ( 1 , 

3 ) - 

PI 3, 

3 ) 


01 132 = SORT ( XT ** 2 4 YT**2 4 ZTv*2) 

ANG = ARC OS ( ABS ( X T*P ( 3, 1) 4 YT*P(3, 2) 4 ZT*P(3, 3))/ 0102) 

DEG = ANG / .0 1745325 3 

WRITE (6, 77) 13, J3, 11, Jl, OES, ANG 

LINES = LINES 4 1 
IF (LINES .LT. PAG SIZ) GC TO 17 
WRITE (6, 75) 

LINES = 1 

17 CALL A V ( P ( l , l), P(l, 2), P(l, 3), P 1 3 * 1), P(3, 2), PC3, 3), 

* P(2, 1), P( 2 , 2), P ( 2, 3), AREA, VOLUME) 

IF (2 * Jl .LE. 11) GC TG 16 

AREA = .5 * AREA 
VOLUME = .5 * VOLUME 

18 SUM OF A = SUM CF A 4 AREA 
SUM OF V = SUM OF V 4 VCLLME 
DO U INDEX l = 2, 4, 2 

XA = P( 1, 1 ) - F( INCEX l, 1 ) 

YA = P( 1, 2) - P( INDEX 1, 2) 

ZA = P( 1, 3) - P( INCEX 1, 3) 

XR = F(3, 1 ) - P< INCEX l, 1 ) 

YB = P(3, 2) - P( INCEX 1, 2 ) 

ZB = P ( 3, 3) - P ( INCEX 1, 3 ) 

0102 = SQRT((XA**2 4 YA**2 4 ZA**2)) * SJR ( ( ( XB**2 4 Y8**2 4 

* ZB**2 ) ) 

ANG = ARCOS((XA*XB 4 YA*YE 4 ZA *Z6> / D102) 

DEG = ANG / .0 1 745325 3 

WRITE (6, 78) II, Jl, LA8(1, INDEX 1), L AB ( 2 , INDEX 1), 13, J3 , 

I OEG, ANG 
LINES = LINES 4 1 


A 




°AG SIZ) Gc .0 19 


IF (LINES .LT 
WRITE ( 6 , 7 5 ) 

LINES = 1 

19 IF ( ! 4 . G T . N ) G J T C 2 1 

CAIL Ay/(P(1, i) t p(i, p (lf 3 > t p (3t 1)f p,, ( j), 

* p <4» l), P(4, 2), p ( 4* 3), AREA, VOLUME) 

IF (2 * J 1 .LE. ID GC TC 20 
AREA = .5 * AREA 
VOLUME = .5 * VOLUME 

20 SUM DF A = SUM OF A 4 AREA 
SUM OF V = SUM OF V * (/GLUME 

21 A l = CtEFIPIl, 25, F(2, 2), P(3, 2), P(l, 3), P(2, 3) 

l ), PI 3, 


ril 

Cl 

A2 

82 

C2 


3 CCEFIPd, 1), F ( 2 


„ - - 1 ) * PU , 

CGEFIPIl, 1), P I 2 , 1), P ( 3 , 1), p(i f 
CCEF( PCI, 2), P ( 4 , 2), P ( 3 , 25, P( 1, 


3), P C 2 , 3) 
2 ) , P (2 , 2 ) 
3), P(4, 3 ) 


3 5 » P14, 3) 


CLEFIPIl, 1), F(4, 1), P ( 3 , 1), p(l, VfF rt , f 31 
CLEF(P(1, l), F{4, 1), P ( 3 , 1), PCI, 2), P(4, 2) 
C1D2 = SQRT ( ( Al**2 4 ei**2 4 Cl**2> * (42**2 4 32**2 
ANG = ARCOS ( -ABS ( A1*A2 4 Bl*B2 4 C1*C2) / C 102 ) 

IF (14 .ST. N) ANG = 2. * (3.1415927 - ANG) 

DEG = ANG / .017453293 

write (6, 79) u, Jl, 13, J 3, DEG, ANS 

LINES = LINES 4 1 

IF (LINES .LT. PAG SIZ) GC TU 23 

WRITE (6, 75) 

LINES = 1 

23 CONTINUE 
WRITE (6, BO) 

LINES = LINES 
IF (LINES .LT 
WRITE (6, 75) 

LINES = 1 

24 CONTINUE 


4 3 

. PAG SIZ) GC TO 24 


A OF Tk I = 2. * SUM C F A 

V OF T S I = 2 . * S JM C F V 

A OF FIG •= F I G < 1, K) ♦ A CF TRI 

7 OF FIG = FIGIl, K) * V CF TRI 

L OF FIG = FIGI2, K) * SUM OF L 

WRITr (6, 81) A OF TRI, V CF TRI, L OF TRI, A CF FIG, 
* L CF FIG 

25 CONTINUE 
GO TJ 1 
END 
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P < U 3), 


♦ P I 3 , 3) ) 
. P I 3 , 3) ) 
, P ( 3 , 2) ) 
» P I 3, 3) ) 

♦ P ( 3 , 3) ) 
. P ( 3 , 2) ) 

♦ C2 **2) ) 


V OF FIG, 


SI 6FTC ROTATE NODECK 

SUBROUTINE ROTATE (N, X2, Y2, Z2. X2R» Y2R» Z2R » XI, fl* ll) 

if <N .EQ. 21 GO TO l 

01 = S CRT ( X 1 **2 ♦ ri«2 ♦ Z 1**2 ) 

HI * XI / 01 

Ui = Yl / 01 

VI * ZI / 01 

H3 * Yl * Z2 - Y2 * il 

U3 * Zi * X2 - XI * 22 

V3 * XI * Y2 - X2 * Yl 

02* SQRT(H3**2 ♦ U3**2 ♦ V3**2) 

H3 * H3 / 02 

U3 * U3 / 02 

V3 * V3 / 02 

H2 * U3 * Zl - Yl ♦ V 3 

U2 * V3 ♦ XI - Zl * H 3 

V2 * H3 * Yl - XI * U3 

03 * SORT ( H2**2 ♦ U2*«2 ♦ V2**2) 

H2 = H2 / 03 

U2 = U2 / 03 

V2 * V2 / 03 

X2R = X2 * Hi ♦ Y2 * Ul ♦ Z2 * VI 

Y2R = X2 * H2 ♦ Y2 * U2 ♦ Z2 * V2 

RETURN 

l X2R = X2 * HI ♦ Y2 * F2 ♦ Z2 * H3 

Y2R * X2 ♦ Ul ♦ Y2 * U2 ♦ Z2 * U3 

Z2R * X2 * VI ♦ Y2 * V2 ♦ Z2 * V3 

RETURN 
ENO 
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S! 


BFTC INTERC NODECK 

SUBROUTINE INTERC <Xi, Yl, ZI 
* X, Y, Z) 


A1 = Y2 - YI 
Bi = XI - X2 
Cl = X l* ( Yl - 
A2 = YA - Y3 
B2 = X3 - XA 
C2 “ X 3 ♦( Y3 - 
A3 = n 
B3 = Zl - Z2 
C3 = Z 1 * ( Y 1 
AA = A 2 
BA = Z3 - ZA 
CA * Z 3* ( Y3 
X - (B2*Cl - 
Y * (A2*Cl - 
Z = I 84*C3 - 
RETURN 
END 


Y 2 ) 4 Y 1 * ( X2 - 


YA) 4 Y 3 ♦ ( XA - 


Y 2) 4 Y 1 * ( Z2 - 


YA) 4 Y3*(ZA - 
B 1*C2 ) / (A2*Bl 
A1*C2) / ( Al*B2 
B3*CA) / ( AA*B3 


t X2, Y2» Z 2 » 


XI ) 


X 3 ) 


Zl) 


Z 3) 

- A1*B2) 

- A2*B1) 

- A3*BA ) 


X3 « Y3 , 


1 3, XA, YA, ZA, 
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M8FTC AV NUDECK 

SUBROUTINE AV (XI, VI , Zl, X2, Y2, Z2, X3, Y3, Z3, AREA, 

XI 2 = XI - X2 

XI 3 = XI - X3 

X23 = X2 - X3 

Y12 = Y1 - Y2 

Y13 = Yl - Y3 

Y23 = Y2 - Y3 

Z12 = Zl - Z2 

Zl 3 = Zl - Z3 

Z23 = Z2 - Z3 

A * SORT ( X 1 2**2 ♦ Y 12 **2 ♦ Z12*-*?) 

8 = SORT ( X 1 3**2 ♦ Y 13**2 ♦ Z13v*2' 

C = SCRT (X23**2 ♦ Y23**2 ♦ 723 **2» 

S = (A ♦ B ♦ C) / 2. 

AREA « SORT (S*(S-AH(S-B)*( S-C) ) 

A * Yl * Z23 - Y2 * Z 1 3 ♦ Y 3 * ZU 

B = XI * Z23 - X2 * Z 1 3 ♦ X3 * Z12 

C * XI ♦ Y23 - X2 * Y 1 3 ♦ X3 * Y12 

0 = XI * ( Y2*Z3 - Y3*Z2 ) - X2 * (Yl*Z3 - Y3*Z1I ♦ 

♦ X3 * (Y1*Z2 - V2*Z 1 ) 

VOLUME = AREA * ABS(C) / SQRT(A**2 ♦ B^2 ♦ C**2I / 3. 

RETURN 

END 


VOLUME) 


Identification Pg lu hedra l Subdivision Identification of Subdivision Method 

JL 


Frequency 6 Icosahedron 


Method 2 


No . of faces 


No. of Uo : of No. of No. of No. of 

vertices sages faces vertices for edges for for total 

“ or PIT PIT f or NIT total sphere total sphere sphere 


Vertex 

Identification 

L 


Vertex 

Identification 


Vertex 

Identification 

L 


1 


1 


63"F(L) = 36' 

'V(G) = 362' 

'E(G) = 1080' 'F 

x coordinate 

, L , 

y coordinate 

z coordinate 

/ 

0.000000 

' 0.850651 7 

; 0. 525731' 

x coordinate 

, L , 

y coordinate 

/ 

z coordinate 

/ 

0.107846 

'0.758171' 

/ 0. 643075' 

x coordinate 

, L 

y coordinate 

/ 

z coordinate 

/ 

0.174499 

'0.66018' 

‘ 0.468576' 


= length 
of element 

'Length' 

angle identification 
angle ft end pt. vertex end pt. 

r-L—r / L 

Axial 0.0 0.0 1 0 1 


Element Identification Unit of length 

^ 1 0 7 '0.21 569298' 


degrees 

± 


radians 


1 '83.808804' ' 1 . 462740' 


angle identification 
ingle a end pt. vertex end pt. 


Face 7 H T 


0 0 10 


■7 / 


degrees 

/ 

71 .646295' 


raaians 

, L 

1.250464 


angle identification 
angle a end pt. vertex end 

J-, 1 L 


aearees 


raj tans 


^ i 21 1 (T 7 '"67.513278' 'l .178329' 
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angle 8 

, — l — 7 

Dihedral 


element identification 
dihedral edge 



degrees 

, — l — 7 

175.363857 


radians 

> L , 

3.060677 


area , volume t length 
for one PIT 

, L 7 

One Face 


Area 


Volume 


Length 

area , volume t length 
for total svhere 

, (— 7 

Complete Structure 


Area 


Volume 


Length 


unit area 

' 0.623048 7 

unit volume 

, L , 

0.206274 

unit length of elements 

' 12 . 489358 ' 


total area 

, L 7 

12.460959 

total volume 

, L 7 

4.125473 

total length 
of elements 

' 216 . 619799 ' 
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frequency 

6 ICOSAHEDRON 

METHOD 2 

V ( L 1 * 28 

El L) = 63 F ( U 

= 36 


362 E(G) * 1080 F< 


0 0 0.000000 


0.850651 0.525731 


i 0 0. 107846 


0.758171 0.643075 


1 1 0.174499 

LENGTH 1 i 

axial 0.0 0.0 

face 1 1 

FACE l 1 

DIHEDRAL 1 1 


0.866018 

1 0 

1 0 

0 0 

2 1 

l 0 


0.468576 

1 1 
1 C 

1 c 


0.21569298 
83 .808804 
71.646295 
67.513278 
175.363857 


2 0 0.212031 


0.639950 0.738585 


2 


0 

.291736 

LENGTH 

2 

AXIAL 

0.0 

FACE 

2 

face 

2 

DIHEDRAL 

2 


0.759807 
1 2 0 

0.0 2 0 

1 1 0 

i 3 1 

1 2 0 


0.581020 

2 1 
2 C 
2 C 


0.2* 341332 
83 .874495 
68.618366 
64.858814 
174.704821 


2 2 0.343074 


0.85 1981 


0.395511 


3 0 0.309017 


0.500000 0.8C9017 


3 


1 0.401810 


LENGTH 3 
AXIAL 0.0 
FACE 3 
FACE 3 
0 I HEDRAL 3 


0.624919 
1 3 0 

0.0 3 0 

1 2 0 

l 4 i 

l 3 0 


0.669347 

3 1 

3 C 
3 0 


0.20910066 
83.998736 
65 .240479 
6 2 • 1C 5942 
174.030886 


3 2 


0.466634 0.729806 


0.459636 


G) * 720 


1.462740 
1.250464 
l. 178329 
3.060677 


1.463886 
1. 157616 
1. 132000 
3.049174 


1. 46 6055 
1.138661 
1.083953 
3.037412 
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LENGTH 3 
AXIAL 0.0 
FACE 3 
FACE 3 
0 IHEORAL 3 


2 3 1 

0.0 3 1 

2 2 ». 

2 4 2 

2 3 l 


3 2 

3 1 

3 1 


0.20977452 

83.979328 

64.955155 

61.944237 

173.958036 


3 


0.500000 


0.809017 


0.3C9017 


4 0 0.395511 


0.343074 0.851981 


0 . 

499636 

0.466634 

LENGTH 

4 1 

4 0 

AXIAL 

0.0 0.0 

4 0 

FACE 

4 1 

3 0 

FACE 

4 l 

5 l 

DIHEDRAL 

4 l 

4 0 


0.729806 

4 1 

4 C 
4 C 


0.20257315 
84.166736 
61 .6C9346 
59.236003 
173.371073 


0 . 

LENGTH 

577350 

4 

2 

AXIAL 

0.0 

0.0 

FACE 

4 

2 

FACE 

4 

2 

DIHEDRAL 

4 

2 


0.57735C 
4 1 

0. 

577350 

4 

1 

4 

2 

3 

1 

4 

1 

5 

4 

2 

l 

4 

1 


0.2 036 1459 
84.150986 
61.316433 
59.027875 
173.262043 


4 


0.624919 


0.669347 0.401810 


4 4 0.639950 


0.738585 0.212031 


5 0 


0.468576 


0.174499 0.866018 


5 1 0.581020 

LENGTH 5 

AX! AL 


FACE 

FACE 

DIHEDRAL 


1 

0.0 0.0 

5 1 

5 1 

5 l 


0.291736 
5 0 

5 0 

4 0 

6 1 

5 0 


0.759807 

5 1 

5 C 
5 C 


0.19408489 

84.431107 

57.869432 

56.243358 

172.752892 


5 2 0.669347 

LENGTH 5 2 

axial 0.0 0.0 


0.401810 
5 l 
5 1 


0.624919 
5 2 


0.19522481 
84 .358302 


1.465716 
1.134380 
1.08 1131 
3.036141 


1.469336 

1.075286 

1.033863 

3.025896 


1.468712 

1.070174 

1.030231 

3.023993 


1.473601 

1.010012 

0.981632 

3.015107 


1.473C28 


•4b- 


FACE 5 2 

FACE 5 2 

DIHEDRAL 5 2 

6 l 

6 2 

5 1 

5 1 

5 1 

57 . 5C242C 
55.9C5183 
172.611000 

1.003607 
0. 975730 
3.0126 32 

3 0*7 29806 

LENGTH 5 3 

axial 0.0 0.0 

face 5 3 

FACE 5 3 

DIHEDRAL 5 3 

0.499636 
5 2 

5 2 

4 2 

6 3 

5 2 

0.466634 

5 3 

5 2 

5 2 

0.19565051 
84 .366042 
57.367121 
55.768106 
172.559191 

1.472814 
1 • CC 1245 
0.97 3686 
3.011726 

’* 0.759807 

0.581020 

0.251736 



> 0.758171 

0.643075 

0.107846 



> 0.525731 

-0.000000 

0.850651 



0.643075 

LENGTH 6 1 
AXIAL 0.0 0.0 
FACE 6 1 
DIHEDRAL 6 1 

0.107846 
6 0 
6 0 

5 0 

6 0 

0.758171 

6 1 
6 C 

0.18426310 
84.713752 
54.176847 
172 .34095C 

1.478534 
0. 94 5564 
3.007917 

! 0.738584 

LENGTH 6 2 

AXIAL 0.0 0.0 

FACE 6 2 

DIHEDRAL 6 2 

0.212031 
6 1 
6 1 

5 1 

6 1 

0.639950 

6 2 
6 l 

0.18426310 

84.713747 

53.51220C 

172.151842 

1.478534 

0.933964 

3.0C4617 

l 0.809017 

LENGTH & 3 

AXIAL 0.0 0.0 

FACE 6 3 

DIHEORAL 6 3 

0.309017 
6 2 
6 2 

5 2 

6 2 

0.500000 

6 3 

6 2 

0.18426311 

84.713752 

53.150170 

172.047482 

1.478534 
0.927645 
3. OC 2795 

0.851901 

0.3955 J l 

0.343074 



0.866018 

0.468576 

0.174499 



0.850651 

0.525731 

-0.000000 




ONE FACE 
AREA 
VOLUME 
LENGTH 


0.623048 

0.206274 

12.489358 


COMPLETE STRUCTURE 

AREA 12.460959 
VOLUME 4.125473 
LENGTH 216.619799 




APPENDIX II 


The plot routine was written for the IBM 7044 
computer utilizing FORTRAN language. The program 
generates data necessary for input tapes to drive a 
Cal Comp plotter model 470. 

The example of the output map is for a six fre- 
quency Icosahedral sphere oriented in three directions 
as shown in Figures 15, 16, and 17. 
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$ ! BFTC PTTOI NODECK 

DIMENSION COORD ( 3 > 3, 32), A(3, 3, 8), 6(3, 3, 8), F(3, 3, 8), 

1 D(3, 3, 8), BUF F ( 4000 ) , ALPHA ( 8 ) , V(3, 3), L AB ( 2, 3), P(3, 3), 

2 SD 1(3, 32), SO 2(3, 32), SD 3(3, 32), EX(768, 3), EXP(768, 3), 

3 C ( 3 , 3) 

COMMON OC ( 6 , 3), X2R2, Y2R2, Z2R2, X2R1, Y2R1, X1R1 
EQUIVALENCE (C00R0(1, 1, I), A), (COORDd, I, 9), B), 

1 (COOROd, l, 17), F), (COORDd, i, 25), D), 

2 (Vd, I), XI), ( V ( 2 , I), Yl), ( V ( 3 , 1), Zl>, 

3 (Vd, 2), X2 ) , t*M2, 2), Y2), (V(3, 2), Z2), 

A (V(l, 3), X 3 ) <V(2, 3), Y3 ) , ( V ( 3 , 3), Z3), 

5 ( LAB ( 1 , 1), d) (LABI2, 1), Ji), (LABd, 25, 12), 

6 (LABd, 2), J2), (LABd, 3), 13), (LABd, 3), J3> 

LOGICAL OUT 

INTEGER RU, TU, RAO, TE OC IC, Q 

COEF ( A I , A2, A3, BI, B2, B3) = A1 * (B2 - B3) - A2 ♦ (BI - B3) ♦ 

1 A3 * (81 - 82) 

DATA A / -.57735027, -.57735027, .57735027, .57735027, 

1 -.57735027, -.57735027, -.57735027, .57735027, -.57735027, 

2 -.57735027, -.57735027, .57735027, .57735027, .57735027, 

3 .57735027, .57735027, -.57735027, -.57735027, -.57735027, 

4 -.57735027, .57735027, -.57735027, .57735027, -.57735027, 

5 .57735027, .5/735027, .57735027, .57735027, .57735027, 

6 .57735027, -.57735027, .57735027, -.57735027, .57735027, 

^ ~ *57735027, — .5 f 735027, 1. ,0.,0.,0.,0.,1.,0.,1.,0.,1.,0«,0.,0., 

8 l.,0.»0.,0. ,— 1. , 1. ,0. ,0. ,0. ,0.,-l. ,0. ,— 1. ,0. , 

*7 1. ,0.,0.,0., — 1. ,0. ,0. ,0., 1./ 

DATA B / — 1., 0.,0. ,0. ,0. , L.,0. , — 1. ,0. , 

1 *'1«,C.,0.,0.,1.,0.,0.,0.,1., 

2 ~~1.,0.,0.,0«,0. ,— l.,0«,l.,0«, 

3 -l.,0.,0.,G.,— l.,0.,0.,0.,-i., 

40. . .8506508 1.. 52573111. .52573111.0. ..85 06 5081. .8 506 5081, 

5. 52 57 3 11 1.0. . 0.,-. 85065081, .525731 11,. 8 5065081, -.52573 11 1,0., 

6.52573111.0. . .8 506508 1.0., -.85065081,. 52573 111, -.52573 111,0., 

7. 8506508 1,-. '85065081, -.5257 31 11,0. ,0., .85065081,. 52 573 111, 

8-. 85065081,. 5257311 1,0., -.525731 11 ,0., .85065C81 / 

DATA F / 0., .85065081, -.52573111, .85065081, .52573111, 0., .5257 
131' 1 ,0 . ,-.85065081 ,0.,-. 8 5065081, -.52 573111*. 52573111,0.,-. 8506508 

2 1.. 8506508 1, -. 525 73 1 11, 0. ,0.,-. 85065081, -.525 73 111, -.8 5065081,-. 52 

35731 1 1.0. , -.525 731 11,0.,-. 85065081,0., .85065081, -.5257 311 1,-. 52573 

4111.0. ,-.85065081 ,-. 85065 08 1, . 52573111, 0.,0.». 85065 081, .525 73 111,- 
5.52573 Id ,0. , .85065031 , .52573 1 1 1,0 .,. 85065081 , 0. ,-.85065081 ,.52573 

61 1 1.. 325731 11.0.. . 85065081, -.52573 111, 0.,. 8 506508 1,0., .85065081,-. 
7525731 ) l, .52573111,0.,-. 85065081, -.5257 311 1,0. ,-. 8506508 1 , 0 850 
86508 l, -.52573 111, -.525 73 1 11,0., -.8 506508 1, .525 73111,0. ,-.85065081/ 

DATA 0 / 0., .8506508i, .52573111, .85065081, . 525 7 31 1 1 , 0 . , 0. , . 850 
16508 1, -. 52573 U 1,0.,. 85065081,. 52 573 11 1,0.,. 85065081, -.5257 31 l l,-. 

285065081. . 52573111 .0. .0. ,-.85065081, .525731 11,0.,-. 85065081 ,-.5257 

3311 1. . 85065061,-. 5257311) , 0 ., 0 .,-. 8 506508 1 ,. 52573 1 1 1 ,-. 85065081 ,-. 
4525 73111,0. ,0. , -. 8 506508 1 , - . 5 257 3 1 1 1 , . 8 506508 1 , . 52 5 73 1 1 1 , 0 . ,.52 573 

51 11.0. .. 85065081. .85065081, -.525 73 1 1 1,0.,. 85065081, .52573 1 11,0.,. 8 

65 06 508 l,-. 52573111, 0. r . 5? 573dl»0. } -. 8 506 5081, -.8 50 65081, .52573111 

7.0. , — .85065081,— .52573111,0. ,-.52573111,0. ,. 8506508 l ,- .8 5065081 , . 5 

82573111. 0. ,-. 52573111,0. ,-.85 06 50 81 ,- . 8506 508 1 , -. 52 573 1 l 1,0. / 

DATA RAD / 1HR / 
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1 

70 


112 

113 


114 

115 


2 


3 

4 


* 


CALL PLOTS (BUFF, 4000) 

CALL PLOT (0., H.5, 23) 

CALL PLOT (0., -11., 21) 

CALL PLOT (6., -.5, 23) 

READ (5, 70) N, TE OC IC, ROT, RU, TILT, TO, ALPHA, METHOD 
FORMAT (212, 6X, 2(F9.0, Al), 8A6, 12) 

IF (N .EQ. 0) STOP 
CALL PLOT (0., 5.5, 23) 

CALL SYMBOL (-5., -4.50, .14, 10HFREQUENCY , 0., 10) 

CALL NUMBER (-0., -0., .14, FLOATIN), 0., -1) 

CALL SYMBOL (-5., -4.75, .14, 9HR0T AT I ON , 0., 9) 

CALL NUMBER (-0., -0., .14, ROT, 0., 3) 

IF (RU .NE. RAD) GO TO 112 

CALL SYMBOL <-0.> -0., .14, 8H RADIANS, 0., 8) 

GO TO 113 

CALL SYMBOL (-0., -0., .14, 8H DEGREES, 0., 8) 

CALL SYMBOL (-5., -5.00, .14, 5HTILT , 0., 5) 

CALL NUMBER (-0., -0., .14, TILT, 0., 3) 

IF (TU .NE. < RAD ) GO TO 114 

CALL SYMBOL '(-0., -0., .14, 8H RADIANS, 0., 8) 

GO TO 115 

CALL SYMBOL (-0., -0., .14, 8H DEGREES, 0., 8) 

CALL SYMBOL (-5., -5.25, .14, ALPHA, 0 ., 48) 

IF (RU .NE. RAD) ROT * ROT * .017453293 
IF (TU .NE. RAD) TILT = TILT * .017453293 
CALL DATE (4., -5., .14) 

CALL FACTOR (5.) 

CALL SYMBOL (-1.2, -1., .05, 3, 0., -1) 

FLOAT N = N 
SIN R * SIN(ROT) 

COS R = COS ( ROT ) 

SIN T = SIN(TILT) 

COS T * COS (TILT) 

MIN = 1 
MAX = 4 

IF (TE OC IC .NE. 8) GO TO 2 
MIN = 5 
MAX = 12 

IF (TE OC IC .NE. 20) GO TO 3 

MIN * 13 

MAX = 32 

DO 4 I * 1, 3 

DO 4 J = 1, 3 

V( J, I ) * COORD ( J , I , MIN) 

LI * 0 

DO 7 L2 * 1, 2 
DO 7 LN = 2, 3 
IF ( L 2 .EQ. LN) GO TO 7 
LI = LI + 1 

CALL RCTATE (1, V(l, LN), V(2, LN), V(3, LN), VII, L2 ) , V ( 2 , 
V ( 3 , L 2 ) ) 

THETA = AT AN2 ( Y2R 1 , X2R1) 


DU 7 L4 = 1, N 
T * L4 

ANG * T * THETA / FLOAT N 


L2), 
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X = COS(ANG) 

Y = SIM ANG) 

A I = -Y 

B 1 = X 

A2 = -Y2R1 

B2 = X2RI - 1. 

C2 = B2 

X = I B 1*C2 ) / ( A2*E 1 - AI*b2) 

Y = (A 1 * C 2 t / (A1*B2 - A2*B1) 
CALL RCTATE (3, X, Y, 0.) 



IFILl - 

2) 45 

, 5, 

45 

SO 

Hit 

L4 ) 


X2R2 


SD 

1(2, 

L4 ) 

= 

Y2R2 


SO 

K3, 

L4> 


Z2R2 


GO 

TO 7 




5 

SD 

2(1, 

L4 ) 

= 

X2R2 


SD 

2(2, 

L4) 


Y2R 2 


SO 

2(3, 

L4) 


Z2R2 


GO 

TO 7 




6 

SO 

3(1, 

L4 ) 

= 

X2R2 


SO 

3(2, 

L4 ) 

= 

Y2R2 


SO 

3(3, 

L4 ) 

r 

Z2R2 

7 

CONTINUE 





Q = 

0 





00 

16 II 

= 1 

t 

N 


00 

16 J1 

= l 

f 

I 1 


12 

= 11 

- 1 




J2 

= J1 

- 1 




13 

= 11 





J3 

= J 2 





OUT 

= .false 

• 


8 

00 

13 LA 

= 1 

» 

3 


LI - LABI It LA) 

LJ = LAB I 2 * LA) 

IJ = LI - LJ 
NJ = N - LJ 

IF (LI .EQ. LJ) GO TO 9 
IF (LJ . EQ. 0) GO TO II 
IF (LI . NE . N) GO TO 12 



P(L A, 

1) 

= 

SO 

3(1, 

NJ) 


P(LA, 

2) 

= 

so 

3(2, 

NJ ) 


P(LA, 

3) 


SD 

3(3, 

NJ) 


GO TO 

13 





9 

IF (LJ 

1 . 

EQ, 

. 0) 

GO 

TO 10 


P ( L A , 

1) 


SO 

1(1, 

LI ) 


P( L A, 

2) 

= 

SO 

H2, 

LI ) 


P ( L A , 

3) 

= 

so 

1(3, 

LI ) 


GO TO 

13 





10 

P ( L A, 

1) 


XI 




P( L A, 

2) 

= 

Y 1 




P ( L A , 

3) 


Z1 




GO TO 

13 





1 1 

P ( L A, 

1) 

= 

SD 

2(1, 

LI ) 


P ( L A , 

2) 


SD 

2(2, 

LI ) 


P(LA, 

3) 

= 

SD 

2(3, 

L I ) 
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c 


c 


c 


GO TO 13 

12 CALL INTERC (SO 1(1, LI), SO 1(2* LI)* 

1 SD 2(2, LI), SO 2(3, LI), SO 3(1, N J ) 

2 SO 1(1, L J ) ? SO 1(2, LJ), SO 1(3, LJ) 

CALL INTERC (SO 3(1, IJ), SO 3(2, IJ), 

1 SD 2(2, IJ), SO 2(3, IJ), SO 3(1, NJ ) 

2 SO 1(1, LJ), SD 1(2, LJ), SO 1(3, LJ) 

CALL INTERC (SO 3(1, IJ), SO 3(2, IJ), 


1 SD 

2(2, 

I J i , SO 2(3, I J ) , SO 

1(1, 

Li, 

2 SD 

2(1, 

LI), SD 2(2, LI), SD 

2(3, 

LI ) 

GO TO 

(1, 

99, 100), METHOU 



99 PILA, 

1) 

* C(l, 1) ♦ C ( 2 , 1) + 

C ( 3 , 

1) 

P ( L A, 

2) 

= C( 1, 2) + C ( 2 , 2) ♦ 

C ( 3 , 

2, 

P ( L A, 

3) 

= C C 1 , 3) ♦ C ( 2 , 3) + 

C ( 3 , 

3) 

GO TO 

13 




100 DO 101 L 

= 1, 3 




SO 1(3, LI), SD2 ( 1 , LI), 

, SO 3(2, NJ), SO 3(3, NJ), 
, C ( 1, 1), C ( L , 2), C ( 1 , 3 ) ) 
SO 3(3, IJ), SO 2 ( 1, IJ) , 

, SO 3(2, NJ), SO 3(3, NJ), 
, C ( 2 , 1), C ( 2 , 2), C ( 2 , 3 ) ) 
SD 3(3, IJ) , SO 2(1, IJ) , 

, SO 1(2, LI), SO 1( 3, LI), 
, C( 3, 1 ),C( 3, 2 ), C ( 3,3 ) ) 


DENOM= SQRT ( C ( L , i)**2 + C(L,2)**2+C(L,3)**2) 
DO 101 M=l,3 


101 C(L,M)=C(L*M)/DENOM 


VERTICES OF WINOOW HAVE BEEN EXPLODED TO SURFACE OF SPHERE. 
DO 102 L= 1 ,2 
DO 102 M= 1,3 


102 C(L»M)=C( L,M)-C(3,P) 

WINDOW HAS BEEN TRANSLATED. 8VERTEX C,I.E. VERTEX 3 IS AT 0,0,0. 
CALL ROTATE (1, C ( 1 » 1), C ( 1 , 2), C(l, 3), C(2, 1), C(2, 2), 

* C(2, 3) ) 

ANG = .5 * ATAN2( Y2R1.X2R1 ) 

XR1 = X2R 1* COS ( ANG ) + Y2R1* SIN(ANG) 

YR1 = Y2R1* COS ( ANG ) - X2R1* SIN(ANG) 

ANG = ATAN2(Y2Ri,XlRl- X2R1) 

ANG = .5 * ANG 

XR2 = (X2R1- X 1 Ri ) * COS( ANG ) - Y2R1* SIN (ANG ) ♦ X1R1 
YR2 = Y2R1* COS ( ANG ) + (X2R1- X1R1)* SIN(ANG) 

51 = YR1 / XR1 

52 = YR2 / (XR2 - X1R1) 

C2 = X1R1* S2 

D9= - Si ♦ S2 

CX = C2 / 09 

CY = (SI * C2 ) / 09 

CALL ROTATE (3, CX, CY, 0.0) 

CENTER OF WINOOW HAS BEEN FOUND. TRANSLATE BACK. . .8 
PiLA, 1) = X2R2* C ( 3 , 1) 

P ( LA, 2) = Y2R2 ♦ C ( 3 , 2) 

P ( L A, 3) = Z2R2* C ( 3, 3) 

13 CONTINUE 

DO 14 INI * 1, 3 

R = SORT ( P ( INI, 1 ) **2 ♦ P ( INI, 2)**2 ♦ P(IN1, 3)**2) 

Q - Q ♦ 1 

DO U IN2 = l, 3 

14 EX (0, IN2) = P ( I N 1 , IN2) / R 
13 = 12 

J2 * J2 ♦ l 

IF (OUT) GO TO 16 

IF ( J2 .GT. 12) GO TO 16 

OUT = .TRUE. 
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GO TO 8 

16 CONTINUE 

WRITE (6, 73) Q 
73 FORMAT ( 5H Q = 14) 

TX1 = XI - X3 
TY1 * Y1 - Y3 
TZ1 = Z1 - Z 3 
TX2 = X2 - X3 
TY2 = Y2 - Y3 
TZ2 = Z2 - Z3 
DO 17 K = 1, Q 
DO 17 K1 * 1, 3 

17 EXP(K, K 1 ) = EX (K, Kl) - V(K1, 3) 

CALL ROTATE (1, TX2, TY2 » TZ2, TX1, TY1, TZi) 

DO 23 K = MIN, MAX 
IF (M .EQ. MIN) GO TO 21 
DO 18 I * 1, 3 
DO 18 J = 1, 3 

18 VI J, I ) = COORD! J, I, M ) 

TX1 = XI - X3 

TY1 = Y1 - Y3 
TZI = Z1 - Z3 
TX2 = X 2 - X 3 
TY2 = Y2 - Y3 
TZ2 = Z2 - Z3 

CALL ROTATE (2, TX2 , TY2 , TZ2, TX1, TY1, TZI) 

DO 19 K = 1, 3 
DO 19 I = 1, 3 

19 P(i', K) = DC ( 1 , K ) *DC ( 4 , I) + DC ( 2 , K)*DC(5, [) + DC(3, K)*DC(6,I) 
DO 20 K = l, Q 

DO 20 L » 1, 3 

20 E X ( K , L ) =EXP(K, 1)*P(L, 1 )+EXP(K, 2 >*P(L ,2 )+EXP (K,3)*P(L, 3) ♦ V(L,3) 

21 00 22 Kl = 1, Q 

YRT = EX ( Kl ,2 ) * COS R - EX(K1,1) * SIN R 

ZRT = EX ( Kl » 3 ) ♦ COS T - EX(K1,1) * COS R * SIN T - EXIK1, :* ) * 

* SIN R * SIN T 
EX ( Kl , 2) = YRT 

22 EX I Kl , 3) * ZRT 

DO 23 Kl « 1, Q, 3 
K2 > Kl ♦ 1 
K3 = Kl ♦ 2 

IF (C0EFIEXIK1, 2), EX(K2 t 2), EX(K3 t 2) f EXIK1, 3), E X ( K 2 , 3), 

* EXIK3, 3)) .LE. 0.) GO TO 23 

CALL PLOT ( EX ( K If 2), EXIK1, 3), 3) 

CALL PLOT IEXIK2, 2), EXIK2, 3), 2) 

CALL PLOT IEXIK3, 2), EXIK3, 3), 2) 

CALL PLOT (EXIK1, 2), EXIK1, 3), 2) 

WRITE (6, 705) (EXINY, 2), EXINY, 3), NY = Kl, K3) 

705 FORMAT ( IX, 6F14.6) 

23 CONTINUE 

call symbol ( ~ i • 2 , —i»* • 0 3 , 3, o«, ~i) 

CALL FACTOR ( 1. ) 

CALL PLOT (8.875, 11., 3) 

CALL PLOT (8.875, 0., 1) 

CALL PLOT (14.875, -.5, -3) 
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HBFTC ROTATE NOOECK 

SUBROUTINE ROTATE IN, X2, Y2, Z2, XL, Yl, Zl) 
COMMON DC 1 6 , 3), X2R2, Y2R2, Z2R2, X2R1, Y2R1, 
IF IN .EQ. 3 > GO TO 1 
Ml = 3*N - 2 
M2 = Ml ♦ 1 

M3 = M£ ♦ 1 

D1 = SORT I Xl**2 ♦ Y 1**2 ♦ Z 1**2 » 

DC ( M 1 , 1) = XI / D1 

DC I Ml , 2) = Yl / D1 

DCIMl, 3) = Zl / D1 

H3 = Yl * Z2 - Y2 * Zl 

U3 = Zl * X2 - XI * 12 

V3-X1*Y2-X2* Yl 

D2 = SCRT(H3**2 + U3**2 ♦ V3**2) 

DC ( M3 , 1) = H3 / 02 
DC I M3 , 2) * U3 / D2 

DC I M3 , 3 ) * V3 / D2 

H2 = DC ( M3, 2 ) * Zl - Yl * DC I M3 , 3) 

U2 = DC (M3, 3) * XI - Zl * DC (M3, 1) 

V2 = DC ( M3 , 1) * Yl - XL * DC(M3, 2) 

03 = SORT ( H2**2 + U2**2 ♦ V2**2 ) 


DC ( M2 , 

1) * H2 / 

03 






DC ( M2 , 

2) = U2 / 

D3 






DC (M2, 

3) = V2 / 

D3 






X1R1 = 

X 1*DC (Ml, 

1 

+ 

Y1*DC (Ml, 

2) 


Z 1*DC (Ml 

X2R 1 = 

X2*DC ( Ml , 

1 ) 

♦ 

Y2*DC (Ml, 

2) 

♦ 

Z2*DC ( M 1 

Y2R1 = 
RETURN 

X2*DC ( M2 , 

1 ) 

+ 

Y2*DC(M2, 

2) 


Z2*DC(M2 

X2R 2 = 

X2*DC ( Ml , 

1 ) 

+ 

Y2*DC (M2, 

1) 

+ 

Z2*DC ( M3 

Y2R2 = 

X2*DC ( Ml , 

2 ) 

+ 

Y2*DC ( M2 , 

2) 


Z2*DC (M3 

Z2R2 - 

X2*DC ( Ml , 

3) 

♦ 

Y2*DC(M2, 

3) 

+ 

Z2*DC (M3 


RETURN 

END 


X1R1 


3) 

3) 

3) 

1 ) 

2) 

3) 
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SIBFTC INTERC NOOECK 

SUBROUTINE INTERC (XI, Yl, ZI, X2, Y2, 


X, 

Y, Z) 





A I = 

Y 2 - Yl 





B I = 

XI - X2 





Cl = 

X 1* ( Y 1 - 

Y 2 ) 

+ 

Y l * ( X2 - 

XI ) 

A2 = 

Y 4 - Y3 





B2 = 

X 3 - X4 





C2 = 

X3*(Y3 - 

Y 4 ) 

+ 

Y3*(X4 - 

X3 ) 

A3 = 

A 1 





B 3 = 

Zl - Z2 





C 3 = 

Z 1 * ( Y 1 - 

Y 2 ) 

+ 

Y1*(Z2 - 

Zl) 

A4 a 

A2 





B4 = 

Z 3 - Z4 





C4 = 

Z3*(Y3 - 

Y4 ) 

♦ 

Y3*(Z4 - 

Z3) 

X * 

( B2*C 1 - 

B1*C2) 

/ ( A2*8 1 

- A 1*B2 ) 

Y = 

( A2*C 1 - 

A 1 *C 2 ) 

/ ( A 1*B2 

- A2*B 1 ) 

Z = 

( B4*C? - 

B 3*C4 ) 

/ ( A4*B3 

- A3*B4 ) 


RETURN 

END 
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— .1 ■ I WfWP * 



Figure 15 


FOLi»Uui I RAMS 2 



FREQUENCY 6 
EO i RT ION 0,000 Lit. GRE 
HOLDOUT FR AXE | TILT 20, 905 DEGREES 
ORIENT R7 T QN C RCE 
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